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Abstract 

Despite their self-reconstruction properties in heterogeneous media, Bessel beams are known to 
degenerate when they are refracted from an isotropic to an anisotropic medium. In this paper, 
we study the converse situation wherein an anisotropic Bessel beam is refracted into an isotropic 
medium. It is shown that these anisotropic Bessel beams also degenerate, leading to confined 
vortical waves that may serve as localized particle trap for acoustical tweezers. The linear nature 
of this degeneration allows the 3D control of this trap position by wavefront correction. Theory 
is confronted to experiments performed in the field of acoustics. A swirling surface acoustic wave 
is synthesized at the surface of a piezoelectric crystal by a MEMS integrated system and radiated 
inside a miniature liquid vessel. The wavefront correction is operated with inverse filter technique. 
This work opens perspectives for contactless on-chip manipulation devices. 
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I. INTRODUCTION 


Bessel beams are helicoidal waves that spiral around a phase singularity (see Fig. [T]) 
where their amplitude cancels resulting in concentric rings aronnd a dark core |il,Q. This 
singnlarity can trap objects and is the cornerstone of aconstic and optical tweezers 
Furthermore, the twisting wavefront transports orbital angular momentum, which can be 
transmitted to absorbing media and exert a torque. This torque can be used to rotate 

fi n n n 

objects [l|, y, llO| or generate controlled vortical flows [ll[. These waves can also be nseful 
for in-depth microscopy since they can self-reconstruct after being damaged by obstacles on 
the propagation axis |3, Q • 


In spite of this healing ability, nnmerous reports 


indicate that refraction between 


isotropic and anisotropic media is fatal to Bessel beams, as crossing the diopter quickly re¬ 
sults in degeneration and fading. In this paper, we show how to predict the fate of anisotropic 
Bessel beams after refraction into an isotropic medinm. Theoretical calculations allows 
the definition of a coherence length, after which the vortex dark core shrinks along the prop¬ 
agation axis, resulting in a dark cavity. As long as a dispersion relation can be set for plane 
waves in both media, our theoretical development allows pre-distorting the wavefront and 
synthesizing confined vortical waves arbitrarily far away from the interface. 

This effect is illustrated experimentally in the field of acoustics. A confined anisotropic 
Bessel beam (called a swirling surface acoustic wave) is synthesized by an integrated 


n 


2Q| at the surface of a piezo- 


MEMS architecture m and inverse filter technique (Ji 
electric anisotropic medium. The refraction of this wave inside a liquid cylinder is measured 
with a needle hydrophone and the resnlt is compared to onr predictions. Then the inverse 
filter method is used for wavefront pre-distortion and displacement of the acoustical trap 
(corresponding to the dark cavity) in the x, y and z directions. 

Vortex confinement and dark cavity 3D displacement are highly desirable featnres since 
they open prospects for dexterous acoustical tweezers. This promising bio-compatible tech¬ 


nology j2]| could manipulate label-free objects jd] with forces up to five order of magnitude 


above optical tweezers. Compared to individual piezo-ceramic elements interdig- 


itated transducers array 


18| has major advantages such as simple integration and low cost 


fabrication with lithography techniques. 
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FIG. 1. Isotropic Bessel beam as defined by equation m with topological order ^ = 1, axial 
wave-vector kz — 1 and radial wave-vector kr — 1. A) Beam cross-section with complex phase and 
amplitude. B) Iso-phase surfaces at W — kzZ — 0 and W — kzZ = vr in red and blue respectively. 
(Color available online) 


II. FROM ANISOTROPIC BESSEL BEAM REFRACTION TO CONFINED VOR¬ 
TICES 


In this section, we will investigate theoretically the refraction of a Bessel beam from an 
anisotropic to an isotropic medium. 

For any physical field /, refraction of waves between two media is ruled by two conditions: 
(i) fulfillment of the wave dispersion relation in each media, (ii) continuity conditions at 
the interface (e.g. continuity of the normal displacement and normal stress in the case of 
acoustic waves). This continuity condition both induces a spatial matching condition on the 
complex wave number and determines the transmitted ft and reflected fr beams amplitude 
resulting from an incident wave /j. We use the superposition principle to resolve / into a 
weighed sum of plane waves. Each of those are defined by a wave-vector associated to a 
complex magnitude describing the field conserved across the interface. According to acoustic 
convention, the locus of the wave-vectors will be referred as the slowness surface, and the 
complex magnitude will be the angular spectrum. 
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A. Dispersion relation condition 


The plane wave decomposition can be simply obtained through an inverse Fourier trans¬ 
form of F, the angular spectrum decomposition of the continuous field /. In Cartesian 
coordinates, we get: 

/ +00 r+OG r+oo 
-oo J —oo J —oo 

where F is the Fourier Transform of /, (x, y, z) are the Cartesian coordinates in the physical 
space and {Kx,Ky,Kz) the Cartesian coordinates in the reciprocal space. Time dependence 
is omitted here for the sake of simplicity. 

Changing variables to cylindrical coordinates gives: k.^ = K.rCos{(f)), Ky = KrSin(0), 
X = rcos(6*) and y = rsin(6*), with (r,9,z) the cylindrical coordinates in the physical 
space and {Kr,(/),Kz) the cylindrical coordinates in the reciprocal space. Straightforward 
trigonometric algebra yields: 

f{r,6,z)= f [ [ F{Kr,(j),Kz)F''^'''^^‘^^^'^~^'^F'^^^KrdKrd(j)dKz (2) 

J — oo J —TV Jo 

For wave fields, the dispersion relation sets a single wave-number k per direction at a 
given frequency: F{K,r,4>, Kz) = h{(f), Kz)S{K,r — kr{(f),Kz)), with kr{4>,Kz) = \JkP‘{(l), Kz) — FI- 
The resulting integral is a general expression for a wave field in an arbitrary material: 

/ —|—^ ~|~oo 

/ h(0, Kz)d(t>dKz (3) 

-TT J —OO 

Now, let ki(0, Kz) be the wave vector of the incident field in the anisotropic medium and 
kt((/), Kz) the transmitted wave-vector. In equation (j31), the angular spectrum h{4>, Kz) is the 
only degree of freedom. For diverse reasons, we can simplify further this equation. First, 
the dispersion relation of the isotropic medium is azimuthal-invariant: 

/ -|-oo r~\~'F 

-OO J -TV 

And second, only the field at the interface is relevant when computing the refraction: 

/ H-oo r-\-7T 

/ ht{(l),Kz)F’'*’''^'^"^''''''''^'^~^ht,r{Hz)d(pdKz ( 5 ) 

-oo J —F 

Similarly to the isotropic medium we only need to evaluate the integral at z = 0 for the 
anisotropic medium: 
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/ “|“Oo 

/ K,)d(t>dK, (6) 

-OO J —7T 

The two integrals (13) and (El) represent the ensemble of waves that could exist in each 
medium depending on the values of hi and hf. They are represented in the reciprocal space 
in figure El hi{(l),Kz) and ht{(f),Kz) are bi-variable functions of complex value on a spherical 
manifold for hf and an arbitrary surface for hi. In order to match the conserved field (e.g. 
normal displacement and normal stress) in both media at a given = 0, we have to balance 
equations (J5l) and db]). This in-plane spatial matching also implies a spectral one in the 
reciprocal space. This means that the projection of both fields on the plane normal to the 
propagation axis must match. The graphical solution then appears as the intersection of 
both media slowness surfaces for all possible 

B. Spatial matching condition 

Besides graphical construction, the refraction of the wave field can be computed with 
the angular spectrum. This framework allows considering each plane wave independently. 
The continuity condition of the physical field / across the interface plane reads fi + fr = ft 
over all the spatial extent of the experiment. Assuming it is much larger than the wave¬ 
length, the same condition applies on the angular spectra hi + hr — ht for any particular 
Kr and (p while keeping a free parameter. We will refer to K,r,z^ for the inci¬ 
dent, reflected and transmitted waves. Using the transmission coefficient t((/),Kz), we can 
express fo(0, Kt^z) = t{<P, K-i,z)hi{(t>, Ki,z)Hki,r{(p, i‘ii,z) “ kt,r{<P, >^t,z)) i with 8 the Dirac operator 
describing the intersection between the slowness surfaces of the isotropic medium and the 
anisotropic one. Hence, integral (Sj) reduces to: 


/ “tOC 

/ Ki^z)hi{(t>, Ki,z)d(t>dKi,z (7) 

-OO J —TT 

The Dirac operator links Ki^z, i^r,z^t,z and enforces a new dispersion relation kzif, Hi^z) = 
K-i,z) = '^\Jkt — kir{(f>, K,i,z)- 25] The locus {ki^r{4>) i^i,z)-: kzif, i^i,z)) describes the slow¬ 
ness surface of any wave that could possibly be transmitted from the anisotropic to the 
isotropic medium. This important result implies that any cross-section of the wave field 
orthogonal to z must fulfill the source-medium dispersion relation. Reversely, given an 
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FIG. 2. Angular spectrum representation of the refraction of a plane wave from an anisotropic to 
an isotropic medium. The incident beam represented by its angular spectrum hi{(j)^i<iz) and wave- 
vector splits into a reflected fleld hr{(l)^tvz) and a transmitted fleld The dispersion 

relation condition appears graphically by forcing all wave-vectors to lie on the slowness surface 
of their respective propagation media. The spatial matching boundary condition forces the wave- 
vectors to share the same propagation direction (p and radial component kr = The 

transmitted (reflected) fleld then happens as the intersection between the isotropic (anisotropic) 
media slowness surface and the set of all wave-vectors with radial component kr. The wave-vector 
k^ is oriented towards the anisotropic medium and cannot exist. Isotropic and anisotropic Bessel 
waves Wi and Wi (dotted blue circles) are the intersection between the slowness surfaces and planes 
normal to the propagation axis. (Color available online) 

arbitrary reference altitude Zq we can synthesize any wave that could be formed on the 
anisotropic medium surface simply by choosing hi{(j)^Ki^z) = 


C. Field refracted by an anisotropic Bessel beam 

In the following, we assume the beam in the anisotropic medium to be either a surface 
or a bulk wave. In the first case, the wave travels along the interface, so hi 0 only for 
= 0. In the second case, we will only consider waves that keep their coherence while 
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propagating along z-axis in the anisotropic medium. For Bessel beams, this latter condition 
amounts to saying that we are able to define some anisotropic Bessel beams that remain 
coherent while propagating along the z-axis. Following this assumption, the incident beam 
angular spectrum must he on the meridian = ki^z- Hence, in both cases hi{(j),Ki^z) = 
hi((l))5{Ki^z — ki,z)- Equation (JTj) becomes: 

f^{r,e,z)= (8) 

J — TT 

Where we used the short-hand notations t(0) = t{(/),ki^z), ^r(0) = and ^^(0) = 

kz{<i>i ki,z)- In our previous work [^, we introduced a definition of anisotropic Bessel beams 
Wi = W['(r, which are coherent solutions in cylindrical coordinates of the wave 

equations in anisotropic media, with Wf the so-called anisotropic swirling wave: 



where I is the topological order of the anisotropic Bessel beam. 

In isotropic media, the integral definition of the Bessel function with a straightforward 
change of variable r] = 7r/2 + (p — 9 yields Wf = leading to the more common isotropic 

Bessel beam IFp 

Wi = ( 10 ) 

Isotropic and anisotropic Bessel beams are defined by their angular spectrum h = (27ri'‘kr{(l>, Kz)). 

In order to create an anisotropic swirling wave \1// of topological order I at an arbitrary al¬ 
titude we set hi{(j)) = /(2m^kr{4>))- Thus, equation ([H|) yields: 


df/(r,6»,z) = 




( 0 )f COs{(l) — 6)-\-ikz{(t>){z — Z' 




( 11 ) 


It is easily inferred from equation (ITTl) that the small azimuthal variations of kz will be 
dramatically magnified by increasing values oi\z — zq\, resulting in beam degeneration. 


III. VORTEX STRENGTHENING, CONFINEMENT AND HOLE FILLING 

For moderately anisotropic media, we can get further insight on this degeneration by 
introducing a simplified slowness surface: ki^r = kf{l + esin(2(/))), with e << 1. The 
computation, detailed in the appendix yields: 
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^i(r, e, z) ~ 0) Jo( 




02 


ez 


z^o 


)+t^e 


.0Ak2z 


+ 00 




n=l 


kf'ez, 

hO ' 


Wf+2jr,ff) + 


-1)”W|-2„(>'.9) 

( 12 ) 


With = kt^ — k^‘^ and the average transmission coefficient over all directions. Equation 
m is made of two terms: a fundamental beam Wf plus a series of degenerated ones W|^ 2 n- 
All these beams are only functions of r, but have a z-dependent factor We 

can note that the same coefficients also appears in frequency modulation signal modeling. 
Indeed, the current phenomenon is highly similar to spatial wavelength modulation since 
the wavelength in anisotropic media is generally direction 0 dependent. For close to 
zero, the zero-order Bessel function predominates, such that the beam is similar to the 
fundamental. For large enough, the Jo(^^) function cancels, so the fundamental mode 
eventually vanishes. In the meantime, the other harmonics and especially the first mode 
grow in amplitude, and reach the fundamental for the first time aX z = K with A solution of 
Jq C\o^ ) = Ji( \o^ ) being the coherence length: 


A ~ 1.43 


We 


( 13 ) 


We can also infer that for even order beams, the degeneration given by equation (1T21) will 
happen through vortex hole filling. With increasing values of n, the terms W°_ 2 n and yVj*^ 2 n 
are magnified by their coefficients J„( \o^ )- Bessel beams of decreasing order have smaller 
dark cores, so for growing n, Wf_ 2 ^ is magnified; it becomes the dominant mode and the 
vortex dark core shrinks. Eventually, provided the vortex is of even topological charge I, a 
large enough value oi n = 1/2 will imprint a Wg wave as local maximum of amplitude on 
the vortex axis, filling the Bessel dark core. This happens at the first local extremum of 


The solution of J^ix) = 0 is often written ^ and is given by Olver 
fi /2 1 = ^(1 + 0.809 X + 0(/“^/^))/2. If we solve d we have: 

Ikl (l + 1.284 X /-2/3) 


^hole — 


2e/cl 


,02 


formula: 


(14) 


According to this equation, the depth of hole of isotropic vortices (e = 0) would be 
infinite. 

To summarize, the following properties are expected when an anisotropic Bessel beam is 
diffracted by an isotropic medium: (i) the vortex can be formed far away from the substrate 
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(ii) the beam degenerates after a distance A which depends on the strength of the anisotropy 
and (iii) even order vortices exhibit hole filling. 

As mentioned in the introduction, these effects are very favorable to acoustical trap¬ 
ping and manipulation. In the following, we test experimentally our predictions and prove 
that confined Bessel beams can be generated from surface acoustic wave integrated devices, 
possibly solving the trade-off between cost and selectivity. 


IV. EXPERIMENTAL SETUP 

The experimental setup, shown in figure ([3D, is similar to the one previously used by 
the authors to synthesize swirling SAW. Surface acoustic waves are generated at the 
surface of a 8 inches 1 mm thick X-cut lithium niobate piezo-electric wafer by an array of 
32 interdigitated transducers powered by a field programmable gate array (FPGA - Lecoeur 
Electronics), whose input signal is designed using two-dimensional inverse filter Q . A small 
PVC cylinder is sealed at the center of the wafer (inner diameter = 14.1 mm, height = 8 mm, 
seal in silicon) and filled with water to study the refraction of the anisotropic swirling wave 
into an isotropic medium. The acoustic field is recorded in water using a 75 /um diameter 
hydrophone (Precision Acoustics LTD) whose position is controlled by a vernier (vertical 
displacement) and 2 micro-stepper motors (horizontal motion). 


Input from FPGA 



FIG. 3. Experimental setup. (Color available online) 
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Each experiment was performed following the same protocol. First, we calibrate our 
32 transducers array by scanning a 2 mm diameter virtual disc located 2 mm above the 
substrate. Then, we compute the optimal input with the inverse filter method to get the 
desired wave field. Neglecting diffraction, there is a unique relation between and ht{(/>). 
So setting the field in a single plane gives hi{(l)), which amounts to choosing in the 

whole fluid cavity. After programming the FPGA, we measure the emitted acoustic field in 
the regions of interest using the needle hydrophone. 


A. Confined acoustical vortex structure 


We synthesized vortices of order I =0, 1 and 2 at Zq = 2 mm above the substrate. Results 
are shown in figure 01 A. 

In figure IHB, we display the numerical integration of equation (fTTI) of the same vor¬ 
tices. The experimental pressure fields are qualitatively and quantitatively similar to the 
predictions for all three topological orders. The experimental field appears slightly blurrier 
than the theoretical one, most likely due to diffraction as the wave meets the three-phase 
interfaces (air/silicone/substrate and silicone/water/substrate) and possibly to the discrete 
number of transducers used in our experiments that limits the synthesis of the prescribed 
wave field. The first bright ring radii scale linearly to the topological order and are respec¬ 
tively 100 and 200 fim for I = 1,2. Despite this lower signal quality, the expected effects 
still appear quite clearly, with the strengthening of the beams, the degeneracy before and 
after the reference plane Zq = 2 mm and the hole filling for the second order vortex. For 
lithium niobate in water at 12 MHz, A ~ 2.2 mm. 

Taking the I = 0 case, we notice that the beam strengthens as it gets further away from 
the substrate. This surprising effect can be explained simply by tracing the origin of the 
beam. When a surface acoustic wave penetrates in the liquid, it starts radiating some power 
Vo along the Rayleigh angle Or ~ 25°. Acoustic waves dissipation in water being very weak, 
these bulk waves reach the vortex axis at z = D/2tan{6R) almost unattenuated. On the 
opposite, the SAW gradually loses power as it propagates under the liquid layer, hence, the 


wave combining at z = 0 will have a power with a ~ 0.20 dB/MHz.mm 27|. 

Following geometrical arguments, the characteristic gain length is /gain = cc;Q;tan(6*/i)/(27r) ~ 
1.1 dB/mm at 12 MHz. Interestingly, some authors proposed that gain media could lead to 
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A. Experimental B. Theoretical 


Amplitude Amplitude 





FIG. 4. Anisotropic Bessel beams propagating in an isotropic fluid. The reference plane is located 
at zo = 2 mm. On this plane, the first bright ring radii scale linearly to the topological order and 
are respectively 100 and 200 /im for / = 1 and I = 2. The beams are expected to change on a 
characteristic distance A ^ 2.2 mm and have their dark core filled after A^oie — 1-8 and 2.8 mm 
for the first and second order vortices respectively. The plotted discs have a 1 mm diameter. A. 
Experimental synthesis of confined acoustical vortices. Topological orders range from 0 to 2. Max 
pressure amplitudes are respectively 0.35, 0.24, and 0.26 MPa. B. Corresponding theoretical waves 
computed from equation f[TT|l . (Color available online) 

tractor beam 

The two other effects are more visible on higher order confined vortices. Looking at 
the phase of the I = 2 vortex, we notice that the phase singularity quickly degenerates for 
increasing \z — Zo\ values and its amplitude bright ring appears squeezed. Similarly, the 
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ring of the first order vortex splits at |z — Zo| > 1 mm. For / = 2, this beam degeneration 
is associated to hole filling, as for increasing \z — Zo\ the beam intensity on the axis grows 
up, filling the vortex hole. This is more easily visualized by taking the cross-section of the 
degenerated vortices, shown in figure [5l 



FIG. 5. Anisotropic second order Bessel beam propagating in an isotropic fluid. The reference 
plane is located at zq = 0 mm. A. Pressure complex amplitude on the vortex axis for various 
altitudes 2 ; normalized by the maximum pressure amplitude of the first bright ring at z = 0. Blue 
circles represent experimental data points and the black solid line the corresponding theoretical 
amplitude. B. Cross-section of the complex pressure amplitude deduced from equation ([TT]) . (Color 
available online) 

The confinement of the vortices suggests to attribute them (x, z) coordinates. In the 
following, we prove experimentally the 3D mobility of these confined acoustical vortices, 
confirming their potential use for acoustic tweezers. 


B. x-y mobility 

As a rule of thumb, larger magnification results in smaller depth of field, so labs on a 
chip are essentially 2D devices. The (x, y) mobility in the plane parallel to the substrate 
is therefore critical. In the present experiment, we used the same calibration to generate 
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vortices at five different {x,y) positions, with the acoustic pressure amplitude and phase 
displayed in figure IHl 



Amplitude 


1 mm 







Phase 


1 mm 



FIG. 6. Synthesis of 2nd order confined acoustical vortices at different positions 2 mm above the 
substrate. Maximum pressure amplitude of each vortex is 0.18 MPa. (Color available online) 

The wave field faithfully reproduces the required displacement over 1 mm (about 15 
wavelengths). This confirms the applicability of confined vortices for 2D manipulation of 
micro-objects. 


C. z mobility 


One of the largest expectation on confined vortices is the vertical manipulation of objects. 
This additional degree of freedom is essential for the manipulation of biological objects like 
cells or bacteria that can strongly interact with walls. This interaction can be desirable to 
achieve cell sorting according to the cells adhesion properties j29l] or detrimental in exper¬ 
iments where walls effects are not the purpose of the study. The z-mobility is also useful 
when studying bubbles or contrast agents properties which are systematically drived close 
to the walls by buoyancy. Finally, it might also allow manipulating multiple stacked labs 
on chips. In this experiment, we kept using the calibration obtained from z = 2 mm, but 
synthesized cross-section of vortices with different reference plane zq ranging from 0 to 4 
mm above the substrate. Results are shown in figure [71 

After the experiment, the data were processed to extract the position of the phase sin¬ 
gularity. The analysis starts by computing the cross-correlation between a confined vortex 
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Target Zo(nnnn) 



Zo = 0 Zo = 1 Zo = 2 Zo = 3 Zo = 4 


FIG. 7. Synthesis of 2nd order confined acoustical vortices at different altitudes zq. A) Maximum 
correlation coefficient. B) Experimental fields. Max pressure amplitudes are respectively 0.22, 
0.24, 0.26 0.30 and 0.33 MPa. (Color available online) 

from equation f fTTl) and the experimental measurements. The maximum of the correlation is 
labeled as the experimental Zq, whereas the one used to design the vortices is tagged target 
Zq. The general good agreement between target and experimental values shown in figure [71 A 
confirms the vertical mobility of confined vortices. The z—mobility is however bounded by 
the spatial extent of the calibration area and the Rayleigh radiation angle. Indeed, waves 
coming from beyond the edge of the calibration area are ill-defined and reach the central 
axis z = Dcai/2 tan(0i^) ~ 2.1 mm away from the calibration plane, bringing some error to 
the wave propagation. 

Visual comparison between vortices of different Zq in figure [TlB agrees well with the 
correlation trend. Here, each column represents a wave with a given target Zq, whereas each 
line is a different measurement plane z. If the vertical motion was ideal, one could expect the 


z (mm) 
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waveforms to be invariant when following the diagonals z — zq = cst. For instance, taking 
the central diagonal {z = zq) we can see that the dark core appears on all scans, whereas 
on the diagonal z — Zq = ±2 mm the hole filling is almost always at the same development 
stage. 

In this experimental section, we successfully synthesized confined vortices at a control¬ 
lable {x,y,z) position. The beams degenerate in excellent agreement with our theory and 
exhibit the researched effect of hole filling. This opens prospects for affordable and dexter¬ 
ous contactless manipulation devices, which will be our future objective. Beyond acoustics, 
we believe depth of hole and vortex confinement could find numerous applications in other 
physical fields involving anisotropic wave propagation. 


V. CONCLUSION 

In this work, we address the general problem of the refraction of an anisotropic Bessel 
beam into an isotropic medium. After getting an explicit expression of the refracted field, 
we derive a coherence length beyond which the beam starts to degenerate. Furthermore, we 
show both theoretically and experimentally that it is possible to anticipate this degeneration 
and form phase singularity arbitrarily far away from the interface. For even topological 
orders, the resulting beam exhibits hole filling. This progressive filling of the vortical wave 
dark core opens prospects for contactless manipulation in both optics and acoustics. In 
optics, similar approach could also allow in depth microscopy in anisotropic media. 
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Appendix: Simplification of the refracted beam for moderately anisotropic materi¬ 
als 


In the following, and without loss of generality, we choose zq = 0. For weak anisotropy, 
we assume the slowness contour of the solid to be ki^r{4>) = ^^(1 + esin(2^)). Using the 
dispersion relation of the isotropic medium + kf,^ = kt^ = we easily get kz{(p) — 
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A:°(l — (A:°/A:°)^esin(20)), with — k^'^. The exponential term e zesin{ 2 (t>)/k° then 

developed according to the Jacobi-Anger expansion, it yields: 




27^^' J-TT feU 

After grouping the terms, we get: 


k^'^ez. 


-A:?^ez, ,1.0^^ ^ .kfez. 


■i',{r,e,z) ~ e^>‘'i‘W';(r,e)M’^)+e'>‘°‘Y.U-^) \WwJr,e) + (-l)"H''L„(r,«) 

n=l 

(A.2) 


With W^i the refracted Bessel beam distorted by the anisotropic transmission coelRcient: 


>V'?(r,0) = ^ r (A.3) 

27rz J—7T 

In practice, we can take for two reasons: (i) errors in complex magnitude on 

the angular spectrum keep constant during the wave propagation, whereas errors in phase 
grow exponentially and (ii) most materials are weakly anisotropic such that the average 
value of t((/)) provides a good approximation of the overall transmission coefficient. This 
yields equation ffT^ . 
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